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Abstract: Prime numbers are natural numbers that can only be divided by one and the 
original number. There is more than one of them. Error-correcting codes used in 
telecommunications are generated using prime numbers. They guarantee automatic 
message correction both during transmission and reception. Algorithms used in 
public-key cryptography are built upon primes. They're also employed in the 
production of pseudorandom numbers. Mathematicians and many other scientific and 
technological communities have always been fascinated by prime numbers. 
Additionally, computer engineers can use it to tackle a wide range of real-world 
problems. The analysis of prime numbers is very important for finding their 
applications in different fields. The statistical analysis of pseudoprimes within a given 
interval is carried out in the presented paper and an algorithm of Python program to 
find the distribution of pseudoprimes has also been generated, which is used to find 
their distribution with different bases within the given intervals. The data analysis 
process made use of graphical depiction. The discovery will surely open up new 
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Introduction 

A composite number x is referred to as a Fermat 
pseudoprime to base a (Wagstaff, 2024) for every 
positive integer a, provided that it divides a*'!—1. The 
Chinese hypothesis is the fallacious assertion that every 
number that passes the Fermat primality test for base 2 is 
prime. 341 is the lowest base-2 Fermat pseudoprime. 
Since it equals 11¢31, it is not a prime, but it does satisfy 
the Fermat primality test for base 2 because it meets 
Fermat's little theorem, which states that 2**°= 1 (mod 
341). When the modifier Fermat is understood, a Fermat 
pseudoprime is commonly referred to as a pseudoprime. 

Sharma et al. (2021) developed a program to locate 
the multi-reverse primes and determine their distribution, 
using a novel definition of primes known as multi-reverse 
primes. Additionally, using multi-reverse primes, a data 
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avenues for future number theory study and applications outside of mathematics. 


encryption and decryption technique has been suggested. 
Sharma et al. (2022a) discovered the correlation between 
innovative primes and defined them, including the digital 
super prime, digital multi-reverse prime, and super twin 
primes. The new primes, known as super primes, were 
developed by Agarwal et al. (2021). A Java application 
has also been created to find different primes within the 
specified range. Pushpa et al. (2021) covered the 
characteristics and outcomes of numbers as well as 
pseudoprime applications 
applications and cryptography. 

The distribution of pseudoprimes is currently an 
important study issue, but it has been challenging to 
analyze their distribution over intervals with long ranges 
because no computer tool has been developed to date. 
This to analyze the distribution of 


various in computer 


study aims 
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pseudoprimes within a specific interval and create a 
Python program to ascertain their distribution. Finding 
the connections between various primes and their uses in 
many scientific fields would be easier with the current 
study's help. 

Literature review 

Somer (1987) established the finite number of d- 
pseudoprimes. The Fermat test was one of the primality 
tests provided (Gradini, 2010). The test's algorithm must 
be coded in Mathematica (6.0 version) in order to be 
performed. The idea of a pseudoprime is introduced by 
discussing the applications of Euler's and Fermat's Little 
Theorems to the tests. Fermat numbers can be used to 
create an endless number of pseudoprimes and super- 
pseudoprimes, as demonstrated by Kfizek et al. (2002). A 
key concept in number theory, Fermat's Little Theorem, 
must be used to define pseudoprimes and _ super 
pseudoprimes. It provides a basic feature of primes and 
serves as the foundation for the majority of primality 
tests. 

He et al.(2022) discussed calculating Fermat's 
pseudoprimes. In 2018, Parthi et al. studied strong 
pseudoprimes, a unique class of pseudoprimes on base b 
that are highly helpful for testing primality. Agarwal et 
al. (2015) generated the function of primes differently 
and created a few distributions as well as their uses in 
cryptography. Agarwal et al. (2023) found a correlation 
between various primes and developed Java software to 
find the correlation model. Agarwal et al. (2021) 
established a new recurrence relation with Fibonacci 
numbers and created a Java program to detect the 
distribution of Fibonacci primes and identify the 
Fibonacci primes inside the specified range. A data 
encryption and decryption technique based on Fibonacci 
primes has also been proposed. Compared to the current 
encryption techniques, the suggested solution offers a 
high level of protection against an unauthorised excess. 

By tabulating all even pseudoprimes up to 101°, 
Pomerance et al. (2023) demonstrated a recent Ordowski 
conjecture, The asymptotic density of the set of numbers 
n that are pseudoprimes to a base that is a suitable divisor 
of n. Li (1996) obtained the upper bound for even 
pseudoprimes in an unpublished study. According to 
Erdos's (1956) conjecture, there are infinite base-a 
pseudoprimes n = r (mod m). If a, r, and m met these 
requirements, Pa;rm(x) of base-a pseudoprimes n = r (mod 
m),n<x is x! as x >, S = {n € N: #D(n) > 0} has 
an asymptotic density, according to the conjecture of 
(Ordowski, 2021); counts up to 10° by A. For sufficiently 
large x, Crm(x) > x!6 '08 loslos 9) if gcd(r, m) = 1. Wright's 
argument for a considerably weaker bound serves as the 
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foundation for this recent result of the first name 
(Pomerance, 1981). In this case, P2:m(x) is unbounded if 
gcd(r, m) = 1. Naturally, this finding from Rotkiewicz 
(1967) is not as strong as the one before it, but it was 
simpler and more than fifty years earlier. 

Sharma et al. (2022b) discovered the relationship 
between numerous primes and the correlation between 
them. In 2018, Agarwal et al. produced new methods for 
number theoretic functions and special numbers that are 
well-defined within a specified range. They also created 
flow charts for these numbers and functions. (Zhang, 
2007) given w’; for t € [13,19]; and given reasons of K2- 
and C3-spsp’s < 10°° to prove wi = w’t< w’’ for t > 12. Wn 
is the smallest strong pseudoprime to the first m prime 
bases (Jiang et al., 2014). For 1 <m <8, the precise value 
of Wm is known. Zhang theorized that wo = wio = Wir = Qui 
after discovering the 19-decimal-digit number Qi, = 
3825123056546413051. Jaeschke (1993) defined wx as 
the lowest strong pseudoprime to all of the first k primes 
chosen as bases. Also described the underlying facts and 
methods used to obtain these conclusions, as well as the 
exact values for ws, Wo, w7, and ws as well as upper 
bounds for wo, wio, and wit. 


Materials and methods 

In the current study, the Python program is used to 
count the Fermat pseudoprimes in a particular interval. 
This information is crucial for figuring out the prime 
pattern and how they relate to one another. Utilizing a 
Python program, the data was collected, and by creating 
line graphs in Microsoft Excel, the distribution of Fermat 
pseudoprimes in various intervals was discovered. 
Distribution of pseudoprimes 

For every base a > 1, an unlimited number of 
pseudoprimes exist. Cipolla demonstrated in 1904 how to 
generate an endless number of pseudoprimes to base a > 
1. In that case, n = AB is a pseudoprime to base ‘a’ and 
composite (Ribenboim, 1996; Hamahata et al., 2007). For 
instance, A = 31, B = 11, and n = 341 are pseudoprimes 
to base 2 if a = 2 and p = 5. Although they are 
uncommon, there are unlimited Carmichael numbers and 
strong pseudoprimes to any base larger than 1. There are 
245 pseudoprimes below one million, 21853 less than 
25109, and three pseudoprimes to base 2 below 1000. 
Below this threshold, 2163 Carmichael numbers exist and 
4842 strong pseudoprimes base 2. 
Program for finding fermat pseudoprimes within the 
given interval 

def ged(p, q): 

while q != 0: 
P4=qg,p%q 
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Figure 2. Distribution of Fermat Pseudoprimes with base 3. 
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Figure 4. Distribution of Fermat Pseudoprimes with base 5. 
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return p 
def is coprime(x, y): 
return gcd(x, y) = 1 
min=int(input("Min. value of interval =")) 
max=int(input("Max. value of interval =")) 
a=int(input("Base =")) 
pce=0 
for x in range(min,max+1): 
res=(a**(x-1))-1 
cp=is_coprime(x,a) 
count=0 
for i in range(2,x): 
if x%i==0: 
count+=1 
if a>1 and count>=1 and cp==True and res%x==0 : 
print(x, "is Fermat Pseudoprime Number") 


pe=pe+l 


No. of Fermat 


Table 1. Distribution of Fermat pseudoprimes within the given interval. 
No. of Fermat 


intervals [1, 10000], [10001, 20000], [20001, 30000], 
[30001, 40000], [40001, 50000], [60001, 70000], [70001, 
80000], [80001, 90000], [90001, 100000] except the 
interval [50001, 60000]. The interval [1, 10000] has 
maximum number of Fermat pseudoprimes with different 
bases. A minimum of five-base fermat pseudoprimes can 
be found at the following intervals: [1, 10000], [10001, 
20000], [60001, 70000], [70001, 80000], [80001, 90000], 
and [90001, 100000]. Furthermore, the distribution of 
Fermat pseudoprimes has no fixed pattern, which 
strengthens its application in the field of cryptography 
because it will be challenging for an intruder to quickly 
identify the pseudoprime utilised in the cryptographic 
process. 


Conclusion 


The current work will aid researchers in their 


No. of Fermat No. of Fermat 


Interval Pseudoprimes Pseudoprimes Pseudoprimes Pseudoprimes 
(Base 2) (Base 3) (Base 4) (Base 5) 

[1,10000] 22 23 47 20 
[10001,20000] 14 11 23 11 
[20001,30000] 4 | 12 8 
[3000 1,40000] 9 5 12 6 
[4000 1,50000] 6 7 17 9 
[5000 1,60000] 3 5 6 8 
[60001,70000] q 4 10 4 
[70001,80000] 3 5 11 2 
[80001,90000] q 7 10 2 
[90001, 100000] 3 4 5 3 


print("Total number of Fermat Pseudoprimes are 
spc) 


=" 


Result and Discussion 
Table 1 represents 

pseudoprimes within the given interval. 

represents the distribution of Fermat Pseudoprimes within 


the distribution of Fermat 
Figure 1 


different intervals with base 2. Figure 2 represents the 
distribution of Fermat Pseudoprimes within different 
intervals with base 3. Figure 3 represents the distribution 
of Fermat Pseudoprimes within different intervals with 
base 4. Figure 4 represents the distribution of Fermat 
Pseudoprimes within different intervals with base 5. 
From the analysis, it is clear that out of Fermat 
pseudoprimes with base 2, 3, 4 and 5, the maximum 
number of Fermat pseudoprimes with base 4 lies in the 


DOL: https://doi.org/10.52756/ijerr.2024.v44spl.010 


investigation of various primes. The suggested Python 
program will assist researchers in determining the 
number of Fermat pseudoprimes that lie within the 
specified intervals with various bases. The finding will 
undoubtedly expand the possibilities for future number 
theory research as well as its applications in other 
domains. Fermat pseudoprimes can also be used in 
cryptography to secure cryptographic 
algorithms. It can also be used to establish a relationship 


generate 


among different primes. 
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